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Abstract. In this lecture I try to explain the basic dynamical processes 
at work in a radiative zone of a rotating star. In particular, the notion 
of baroclinicity is thoroughly discussed. Attention is specially directed 
to the case of circulations and the key role of angular momentum con- 
' servation is stressed. The specific part played by viscosity is also ex- 

, plained. The old approach of Eddington and Sweet is reviewed and 

QQ ■ criticized in the light of the seminal papers of Busse 1981 and Zahn 

1992. Other examples taken in the recent literature are also presented; 
\Q . finally, I summarize the important points. 

o 

i ' 1 Introduction 

o ■ 

"o^ I It is well known since the seminal work of Ivon ZeineJ l|l924l) that radiative zones 
d • of a rotating star cannot be at rest; in other words no hydrostatic solution exists 
| in any frame: some shearing flow must occur. The importance of such flows 
k~j ' comes horn their relation with the abundances of elements at the surface of stars. 
rN Indeed, the interpretation of these observations requires the understanding of the 
■ close relation between rotation and mixing, namely the transport of elements in 
' radiative zones, which, without rotation would be a pure diffusive process. 

As any fluid dynamical process, flows in rotating radiative zones are governed 
by PDE and therefore difficult to evaluate in a simple way. This is the reason 
why much misleading work has been published on the subject, especially under 
the title of the Eddington-Sweet circulation. In this lecture I would like to present 
as clearly as possible the dynamical processes at work in a rotating radiative zone 
and explain why the idea of Eddington-Sweet circulation is a misunderstanding of 
the problem. 
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2 Dynamical ingredients 

2.1 Equations of motion 

Unlike many problems in fluid mechanics where flows are produced by a given force 
field (just think for instance to the case of a river whose water flows down thanks 
to gravity), flows in rotating radiative zones are the result of a mismatch between 
two force fields, namely the pressure gradient and the gravity force. This mismatch 
gives birth to a torque field, the baroclinic torque, which locally produces vorticity 
and thus fluid motion. Before discussing this matter, it is however necessary to 
have a look to the equations of motion and especially to the momentum equation 
in a rotating frame; this equation reads: 

r | + 2£LW+ f£s ) = V P + />ij 
~~i~~ n / m 

where we discarded viscosity. Terms I, II, III are specific of this problem; they are 
respectively for: 

• The Coriolis acceleration 2ft A pv which insures conservation of angular mo- 
mentum. 

• The centrifugal acceleration —fl 2 s which is at the origin of baroclinicity. 

• The gravity force which transform a density fluctuation into motion. 



2.2 The wave zoo 

Before diving into the heart of the subject, it is useful to recall the waves able to 
propagate in a rotating radiative zone. There are three sorts of waves of that kind: 
the acoustic, the gravity and the inertial waves. 



2.2.1 The acoustic waves 

In a uniform medium acoustic waves are described by 

' Pof = -W 

< ^L+ Po dWv = (2.1) 
I at u 

namely momentum, mass and entropy conservation. Note that for such waves dtp' 
is an essential term. Equations have been linearized and 0-index terms refer to the 
background state. The dispersion relation of these waves, namely u — c s k, where 
c s is the speed of sound, shows that acoustic frequencies have no upper limit. But 
they have a lower limit given by a k min coming from the finite size of the container 
(the star !). 
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2.2.2 The gravity waves 

While acoustic waves exist because of the pressure field, the restoring force of 
gravity waves is the buoyancy. Compressibility is unnecessary and therefore den- 
sity can be kept constant except when it yields buoyancy (this is the Boussinesq 
approximation in a very short description). In its simplest form, the system of 
equation describing these waves is 

Pof = -Vp' + p'g 

f™ = °, (2.2) 
p'/po = -aT' 

where a is the dilation coefficient. Note that temperature, which could be replaced 
by any other scalar field (like concentration for instance), needs to be non- uniform. 

Unlike acoustic waves, gravity waves suffer an anisotropy imposed by the di- 
rection of gravity and temperature gradient. This results in a dispersion relation 
of the form: 



LU = N 



fc2 _l_ 



where gravity and VT are along e z ; N = ( — ag^-^j is the Brunt- Vaisala 
frequency. This relation shows that the wave frequency is bounded by N and 
therefore gravity waves occupy the low-frequency band. 

2.2.3 The inertial waves 

Lastly, waves induced by the Coriolis force. Inertial waves are the solution of 



p (f + 2QAtf) =-Vp' 
divv = 



(2.3) 



and thus contain only a velocity and a pressure perturbation. Like the gravity 
wave they behave anisotropically. The direction of the rotation axis gives the 
preferred direction. Plane waves obey the following dispersion relation: 

k 

which shows that these wa yes are also in the low-frequency r ange, bounded by 2f2 
the Coriolis frequency (see lOreenspanl ll969L iRieutordl Il997l) . 

2.2.4 Some more waves 



One often meets, especially in the geophysical literature, other waves like Rossby, 
baroclinic... waves. These are subclasses of the foregoing (gravity and inertial) 
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waves meeting some additional constrains or evolve on some specific flows. For 
instance, Rossby waves are inertial waves which do not (or little) depend on the 
coordinate parallel to the rotation axis, while baroclinic waves combine gravity 
and coriolis forces and surf on shear flows. 

2.2.5 Further remarks 

As inertial waves and gravity waves occupy the low-frequency range, they can 
perturb each other very strongly. The common situation in stars is that N > 2f2. 
This means that gravity modes of frequency near 2f2 or lower, are strongly per- 
turbed. In that case, gravity and inert ial modes are undisti nguishable and one 
speaks about gravito-inertial waves fsee iDintrans et al"lll999|K 

When dealing with rotation we did not mention the centrifugal acceleration; 
this is because this field does not generate any specific wave. It modifies an existant 
gravity field and therefore just enters the buoyancy force. 

3 Baroclinicity 

3.1 Introduction to baroclinicity 

Baroclinicity means the inclination of isobars with respect to "horizontality" usu- 
ally given by equipotentials. More generally, if isotherms, isodensities, isobars and 
equipotentials are not identical, the fluid is said in a baroclinic state while if iso- 
densities and isobars (and thus isotherms) are identical the fluid is said barotropic. 
The importance of this distinction comes from the existence or inexistence of hy- 
drostatic solutions. Let us examine the necessary condition for a fluid to be in 
hydrostatic equilibrium: 



—VP — pV$ = Mechanical equilibrium 
div(xVT) = Thermal equilibrium 
f(P,p,T) = Equation of state 

V$ is either prescribed or obtained from A$ = AirGp — 2fl 2 . In order that 
mechanical equilibrium be possible, it is necessary that: 

Rot (-VP) = ^> Vp A VP = 



i.e. that isobars and isodensities are identical surfaces. Now we may differentiate 
the equation of state, 

dP dP+ dp- dp+ df dT = ° 

Considering that mechanical equilibrium is realized, on an isobar (which is an 
isodensity), we have 
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Fig. 1. Inclination of an isobar, an isotherm and an equipotential in a baroclinic set-up. 



which shows that either the isobar is an isotherm, dT = 0, or ^ — i.e. the 
fluid is barotropic and there is a relation between pressure and density, P = P(p). 
This latter case is rather specific and usually found because either the fluid is 
considered to be isentropic or isothermal. In general, ^ ^ and isotherms, 
which are determined by the equation of thermal equilibrium, are not identical 
with isobars. This general case is illustrated on figure 1. 



3.2 Examples 

To illustrate the notion of baroclinicity we shall consider two classical examples 
which can be solved explicitely. 



3.2.1 Double glazing window 

The first example comes from the classical problem of natural convection which 
arises within a fluid imprisoned between two vertical plates at different tempera- 
tures; such a situation is the one prevailing inside isolating windows with double 
glazing where air is trapped between two glass panes, a cold and a warm one. 

Since temperature fluctuations are small, we assume the Boussinesq approxi- 
mation and neglect all density fluctuations except those generating the buoyancy. 
Such a small amplitude motion is governed by: 

-V5P + Spg + pAv = 
v ■ VT eq = nAST 
< divv = (3.1) 




= -aST 



Recall that in this problem VT eg = since at equilibrium temperature needs to 
be constant. 

Let us assume that the fluid is between two plates, infinite in the y and z 
directions separated by d in horizontal x-direction, with the warm plate at T = T c 
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(resp. cold at T = Tf) in x = —d/2 (resp. in x = d/2). Thermal equilibrium leads 
to 

5T(x)=-(T c -T f )x/d 

a{T c - Tf)p 
5p= x 

and thus the velocity field verifies: 

fj,Av — VSP = a{Tc -p )pa9 xe z , diyv = 
d 

in this problem the baroclinic torque is balance by the viscous torque: 

hAlS = — — — , IJ e x A e z , w = Vxw 
d 

The solution of this equation is easily derived and may be expressed as: 

vM = a[ \^ P09 <<?-^) (3-2) 

where we used no-slip boundary conditions on the plates together with a zero- 
vertical mass flux to determine the three constants of integration. The flow thus 
obtained is shown in figure [5] 

This is a vortical shear flow where vorticity is generated by the torque and 
diffused to the walls by viscosity. 



3.2.2 Thermal wind 

The second example is more closely related to our astrophysical concern. It is 
issued from geophysics and describes the baroclinic situation that occurs in the 
Earth's atmosphere where solar heating imposes a latitudinal temperature gradi- 
ent; poles which are cooler than the equator may be compared to the cold pane of 
the foregoing example. If the Earth were not rotating the baroclinic torque would 
induce, like above, a circulation between poles and equator. However, Earth ro- 
tates and therefore air cannot move easily from pole to equator because of angular 
momentum conservation. In this case the baroclinic torque induces an azimuthal 
flow; the mechanism is just similar to the one of a precessing spinner: the torque 
induced by the gravitationnal field on the spinner induces an azimuthal motion, 
namely the precession. 

Let's have a look to the equation of the flow using the Boussinesq approxima- 
tion; the steady linearized inviscid momentum equation reads: 

2VlAv = -—WSP-aSTg 
Po 

where we neglected viscous force and non-linear terms. Here the baroclinic torque 
is balanced by the torque issued from Coriolis acceleration, namely 
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Fig. 2. Flow between two vertical plates at different temperatures. This picture illus- 
trates the shear flow generated by the baroclinic torque. 



Rot(2ft A v) = Rot(ag5Te z ) 

More explicitely, assuming that 5T = ST(x) (x being the North-South direction), 
we find 



8Vy 

dz 



oc 6T'(x) 



This equation shows the appearance of an azimuthal v y flow and an axial gradient 
of angular momentum. 

Anticipating the next section, we see that the first effect of a (baroclinic) torque 
is, as in the precessing spinner, to change the angular momentum distribution and 
not to induce a meridian circulation. 



4 Rotating radiative zones - 1924 - 2004 : the evolution of ideas in 
eighty years 

4.1 von Zeipel theorem 1924 

In 1924 von Zeipel addressed the question of the equilibrium of a rotating radiative 
zone. Basically, he considered the point of section A3. II in the stellar context. Let 
us remind us his result. If we assume that a radiative rotating envelope is in 
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hydrostatic equilibrium, then, as we showed in 13.11 isobars, isochores, isotherms 
are identical and therefore all thermodynamical quantities should depend only on, 
say, the total potential $ (centrifugal and gravitational); if the thermal equilibrium 
is realized 

div(xVT) + pe = or div(x($)VT($)) + pe = 



<S=> xT'A$ + (xT')'(V$) 2 +p£ = 

where e is the rate of energy generation per unit mass. Using Poisson equation for 
the effective potential, 



y($)r'($)(47rGp($) - 2Q 2 ) + (x($)T' '($))' ^ (V$) 2 +pe = 

Cte Cte not Cte 

where V$ is just the effective gravity which is not constant on an equipotential 
(Spoio 7^ <?cq)- With such an equation, von Zeipel correctly concluded that an 
equilibrium solution is possible only if — and if 

on * 



2nGp, 

Apparently, von Zeipel was pleased with such a result. But lEddinetonl <|l925|) 
was not, being skeptical that energy generation, presumably of microscopic origin, 
could follow such a law. He suggested 1 that hydrostatic equilibrium should be 
abandoned and that a meridian circulation would result from the thermal imbal- 
ance, putting forward the picture that hot poles and cool equator (or vice versa) 
would lead to some global convection with rising material at poles and sinking one 
at equator (or vice versa). 

Later, the solution of von Zeipel was qualified as paradoxical by lOpild lll95lh 
because of being not physically sound. 1 suppose that it is since that time that 
von Zeipel result is also quoted as "von Zeipel paradox". Presently however, von 
Zeipel paradox rather refe rs to the inexistence of a static radiativ e equilibrium in 
a uniformly rotating star l|Bussel I1982L lHansen fc Kawalerl Il994|) . 

4.2 The Eddington-Sweet answer to von Zeipel paradox 

The ideas of Eddington and Vogt have been transformed into a model by Sweet in 
1950. The reasoning is the following: since div(xVT) ^ and since we are looking 
for a steady solution, the velocity field which appears must be such that: 

pc v (v- VT) = div( X VT) 



1 Independently, Vogt 1925 had the same idea. 
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Assuming the thermal imbalance, div(xVT), is given, the temperature field be- 
ing almost spherically symmetric, this equation gives an expression of the radial 
velocity: 

dT 

pc v v r — ~ div(xVT) 
Or 

while vg is obtained from mass conservation: 

i. / l<9/2 \ 1 9 , . „ , 

dxv(pv) = —^-{r pv r ) H — — (s m 6pv e ) 

r z Or r sin 06 

The problem is that such a velocity field is not produced by a force field and 
therefore the conservation of angular momentum is not insured since it is not a 
solution of the momentum equation. 

Such a difficulty was already pointed out bv lRanderd l|l94lj) who underlined the 
key role of vis cosity and ques t ioned t he existence of meridian curre nts in an inviscid 
star. Then, ISchwarzschildl l|l947|) . iBaker fc Kippenhahn] l|l959f) and iRoxburghl 
( 1964| showed more and more explicitly that some differential rotation permits 
the existence of (inviscid) rotating radiative zone without meridian circulation 2 . 
Th ese authors actually pointed o ut the existence of the thermal wind solutions. 

iKippenhahn fc WeieertJ |l990) in Stellar structure and evolution note the prob- 
lem of angular momentum conservation and the existence of solutions without 
meridian circulation but do not show the link between the two properties. 

In fact , the connection betwee n differential rotat ion, viscosity and meridian 
circulation were explained by Busse <)Busselll98lUl98^ and the first s elf-consisten t 
solutions (up to a turbulence model) were given in a seminal paper bv lZahnl f 1992 V 

We need now to review these two landmarks of the subject. 



4.3 Busse 1981 : "Do Eddington-Sweet circulation exist?" 

The problem faced by von Zeipel was presented in terms of a thermal imbalance 
which, through buoyancy (hot poles and cold equator), would drive a meridian 
circulation. However, such a steady motion is impossible in a rotating fluid be- 
cause the angular momentum is not locally conserved: the flux of angular momen- 
tum induced by the meridian circulation need to be compensated by another flux 
(through magnetic or viscous forces for instance). 

Such steady motion are thus impeded mechanically if we ignore all other forces 
than buoyancy. Another possibility is of course to relax the steadiness hypothesis; 
however, in such a case conservation of angular momentum appears as a restoring 
force (namely Coriolis force) which will transform any initial perturbation into a 
(inertial) wave motion whose transport properties are another story 3 . 



2 Actually, Milne, quoted by EddinetmJ ll92Bl) . also proposed such a possibility. 
3 Wavcs, although being an essentially periodic motion, can transport momentum or chemicals 
through their mean motion when they are of finite amplitude. 
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Fig. 3. Thermal wind in the cartesian analog of rotating radiative envelope. 

However, if the situation is stuck mechanically, it can evolve thermally. This is 
what occurs indeed. In this way the fluid may move... This remark underlines the 
real nature of the problem as it has be enounced: this is an initial value problem; 
then, if the time scale of the transient is short enough the star reaches a steady 
state or, if it is too long, initial conditions should be of dramatic in fluence . 

To illustrate the foregoing discussion, we shall follow the work of lBussd l)l98ll) 
and his example in cartesian geometry. With such a set-up, one can consider 
arbitrary initial conditions and compute the flow and its time scales in every 
detail. 

Busse considered a fluid with very small Prandtl number (as expected in a 
radiative stellar envelope) so that viscosity can be neglected in the first instance. 
Thus, the only source of (long term) evolution is thermal diffusion which damps 
out the transients to the steady state. According to our foregoing remarks which 
pointed the impossibility of steady meridian circulation without viscosity or mag- 
netic field, we expect that such meridian flow occurs only during the transient. 

The configuration is shown in figure [31 to keep things as simple as possible 
Busse considers an axisymmetric flow whose cartesian analog is: 



where A = — a + iuj; considering the flow is of small amplitude all nonlinear terms 
are neglected; the equation are further simplified by using the Boussinesq approx- 
imation. Although the system seems to be much simpler than real stars, the heart 
of the problem, as discovered by von Zeipel, is conserved. Baroclinicity, as it is 
called, is still there. Casting the equations of motions into a single scalar equation 
yields: 



which can be solved by Fourier expansions in the horizontal direction. Solutions are 
split into different sets of damped modes. The main ones are the gravito-inertial 
modes: 




Xt 




Michel Ricutord: On the dynamics of radiative zones in rotating stars 11 



"J/ = cos(mkx) sin (nir(z/d + 1/2)) 

V = — ; coa(mkx) cos (mr(z/d + 1/2)) 

= sin(mfca;) sin {mr{z/d + 1/2)) 

.2 _ 



4fi 2 n 2 + N 2 (mkd/ir) 



a = nk 2 m 2 „ 

2w^ 



n 2 + (fcd/7r) 2 
TV 2 



and the baroclinic modes: 



= sin(m/ca:) sin (nn(z/d + 1/2)) 
mkpagd 

V = cos{mkx) cos {nir(z/d + 1/2)) 

2liK7T 
* = O^/d 2 ) 

w 2 = 

k n 2 7r 2 + m 2 k 2 



d 2 1 + N 2 /m 2 {mk/nir) 2 

These sets of modes should be completed by the geostrophic modes and pure 
thermal modes. 

The fact that the Fourier basis is complete implies the completeness of the set 
of modes and therefore the possibility of decomposing any initial condition into a 
linear combination of the modes. Thus the time scale of the transient is controled 
by the least-damped modes. To obtain the corresponding spherical analog of this 
model we should make the correspondance 

d~R, k~2/R 
Two damping rates associated with the two (main) sets of modes appear 

4K (7T 2 +4)7T 2 k 
(TGI ~ "^2, Tbaro 1 ^ l 4 ' 1 ) 

where rj — < 1. Thus we find two time scales: 

• Tkh = -7- = The thermal diffusion time or Kelvin- Hclmoltz time scale 

• T ED = T KH /r] = The Eddington-Sweet time scale. 

For rapid rotators 77 ~ 1 and all transient are damped on the Kelvin-Helmoltz 
time, while for slow rotators rj <C 1 baroclinic modes decay very slowly and initial 
conditions play a crucial role. 
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What about meridian circulation ? In this model it appears through We 
see that there is some amplitude associated with gravito-inertial modes, but this is 
oscillatory motion which disappears after some diffusion time. Rather, the decay 
of baroclinic modes, which are characterized by no oscillation may be considered 
closer to what was imagined by Eddington, Vogt and Sweet. However, we see that 
either the star is rapidly rotating and such flows disappears on a Kelvin-Helmoltz 
time or it is slowly rotating and it keeps a strong memory of the initial conditions. 

Finally, to make the discussion complete, Busse considered the Ekman layers 
which, through Ekman (viscous) pumping, can induce some circulation, but con- 
sidering the small value of viscosity, such flows turn out to be even weaker. His 
paper ends with a discussion of the stability of the steady flow (the thermal wind) , 
a stability which depends on the amplitude of course. The scenario proposed by 
Busse may be illustrated as in figure 0] 

Back to the title of Busse's paper "Do Eddington-Sweet circulations e xist?" we 
clearly see that the answer is negative. In a following paper (|BusseU l982'l aimed at 
popularizing this idea, Busse hammers it in: "... there are no special limits in which 
the Eddington-Sweet theory provides a correct solution of the basic equations" . This 
strong conclusion comes after forty years (since the work of Randers) during which 
the main stream literature has been forgetting basic physics principles ! 

In fact, the work of Sweet correctly derived the now well-known 'Eddington- 
Sweet' tim e scale , which is the time it takes for the longest transient to be damped 
out. After ISweeH (,1950), this time scale is also understood as the time it takes for 
a fluid element to cross a distance like the radius of the star; we shall see below 
( H4.6I) that this can be true only in an unrealistic case. 

The work of Busse however considers an idealized situation where the transport 
coefficient are of molecular or radiative origin. Hence, the Prandtl number is 
very small. Actually, even a radiative zone may be the seat of some turbulence 
generated by (non-axisymmetric) baroclinic instabilities. Indeed, the thermal wind 
is basically a shear flow prone to shear instabilities. We thus expect that some 
turbulence takes place and enhances strongly the diffusion of momentum in stellar 
plasma. 

But, as we pointed out, the presence of viscosity (or any local diffusive trans- 
port) changes much the situation since it permits the existence of a st eady meridia n 
circulation. We shall now examine this question through the work of lZahnl l|l992fl . 



4.4 Zahn 1992 

The work of lZahnl l)l992|) is basically the implementation of Busse's approach in 
the framework of spherically symmetric stars but taking into account turbulent 
viscosity. 

The main idea is to assume that, thanks to the stable radial stratification of 
a radiative zone, turbulence is essentially horizontal. Hence, turbulent diffusion is 
strongly anisotropic, being most vigorous in the horizontal direction, implying a 
differential rotation developed preferentially in the radial direction. 

The model concentrates on three equations: 
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Rotation + Stratification 



Baroclinic temperature field 
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Baroclinic torque (mechanical imbalance) 



Thermal wind (differential rotation) 



Shear 
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Meridional circulation 



Fig. 4. Busse scenario; see comments in fig. 



• The advection-diffusion balance of the angular momentum flux: 

1 d ^ 4 + 1 d / 4 dn\ _ Q 
5r 2 dr r 2 dr \ dr j 

where v r — U(r)P2(cos 9); this form of the velocity field retains only the first 
term in a Legendre polynomial expansion. 



• The balance of torques, source of the thermal wind, 

= Vft 2 A se s 



VpAVP ^ o2 



P 2 



• The advection-diffusion balance of the thermal flux 
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div(xVT) + pe = pc v v ■ VT (4.2) 

These equations are of course completed by the equation of mass conservation. 

Zahn then imposes an important massage to the equation of energy (|4.2ll in 
order to obtain the expression of the thermal imbalance div(xVT) + pe in terms 
of its barotropic and baroclinic components. Each of these components are, like 
the velocity field, expanded onto Legendre polynomials and only the P2(cos6)- 
component is retained; the thermal imbalance is thus written: 

div(xVT) + pe = P^{En + £ M )P 2 (cos0) 

where Eq and are non-dimensional functions of rotation and compositional 
gradient respectively. 

In section §3.5 of the paper Zahn discusses the case of the meridian velocity 
since its expression is easily derived fro m 114.21) . The re, he s hows h ow to recove r 
the expressio ns formely derived by l»attonl 111941 . ISweetJ l|l95^ . lOnikl lll95lh . 
iMestell l|l966|) when the rotation is assumed uniform. Of course such a rotation 
cannot exist by itself and need to be imposed by some trick (a magnetic field for 
instance although in such a case even the c irculation mig ht be suppressed, see a 
discussion of the case of magnetic stars in iMesteH ((^003)). However, as shown 
by Busse, if the trick is removed, these expression of the meridian circulation are 
pure non-sense. Indeed, a state of rigid rotation is never reach and terms coming 
from the ^-gradients will just cancel the Eq term ! In a steady state, radial 
velocity vanishes if there is no viscosity. The casual reader should therefore not 
overlook section §4.1 of Zahn's paper where the dependence between circulation, 
viscosity and diffe rential rotation i s give n explicitly (see equation 4.6). This point 
is emphasized in lMaeder k, Zahnl 1)1998): "... not only but also Eq evolves to 
stop the circulation in a star which does not lose angular momentum" . 

4.5 A further example 

In order to further illustrate the relations between circulation, viscosity and dif- 
ferential rotation, we shall focus on another simple model, in the style of Busse's 
but which retains viscosity and spherical geometry. 

The set-up is that of a rotating, self-gravitating, full fluid sphere of nearly 
constant density so as to be allowed to use the Boussinesq approximation. The 
boundaries of the domain remain spherical whatever the rotation speed. 

Within such a set-up, fluid motion, in the rotating frame, is governed by the 
following equations: 

2Q A v + v- W = --V8P + — (g + Q 2 se s ) + isAv 
P P 

v ■ V5T = kAST + Q 
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Meridional circulation 



viscosity 



Differential rotation 
Shear 



Baroclinic temperature field 
Baroclinic torque (mechanical imbalance) 



Fig. 5. Zahn's approach. This is the same as Busse's one but put in another ordre. Each 
arrow symbolises an equation; therefore their orientation has no meaning and could 
be reversed! however, it is often seen in the literature that thermal imbalance implies 
meridian circulation. Such a statement is pure non-sense if isolated from the others; wc 
have seen that it would be much better to say that "differential rotation implies, through 
viscosity, meridian circulation". In fact, this diagram just shows that we face a system 
of coupled partial differential equations. 



div v = 

5 -l = -aST 
P 

where v is the velocity field, Q — Qe z the angular rotation of the frame, g the local 
gravity, v the kinematic viscosity, k the thermal diffusivity, a the thermal dilation 
coefficient. 

Using R, the radius of the sphere, as the length scale, (20)" 1 as the time scale, 
as the temperature scale, the non-dimensional form of the equation is: 
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e z Au + u ■ Vu — — Vp — Ra0(— re r + ese s ) + EAu 

u-ve = e t aq + q 

div u = 
with the non-dimensional numbers 

aT.«g Q 2 R 
Ra= 4^' £= — 

which control the driving and the other numbers 

v „ K 



(4.3) 



E 



2QR 2 ' 



En 



2ilR 2 ' 



which control the diffusion. 

The system (|4.3|) is simplified by elimination of the barotropic component of 
the background which is spherically symmetric and in hydrostatic equilibrium; we 
thus write: 



u = eui, 6 = 8o(r) + e6i, p = po(r) + epi 
Furthermore, if we introduce the Brunt- Vaisala, frequency 

Nj- = a—y—g = 4il 2 Rar0g(r) 
dr 

and set HarQ' = rn 2 {r), the equation of motion may be cast into: 



Rot e z A u + eu ■ Vu — (re r — £se,,)8i — EAu 

n 2 u r + eu ■ V6i = E T A6i 
div u = 



-rn sin U cos ( 



(4.4) 



The equation of vorticity clearly shows that the motion is driven by the baroclinic 
torque, —rn 2 sin 9 cos 0e v issued from the set-up. We are therefore facing the prob- 
lem of a steady forced flow. We may further simplify the problem by considering 
small e, i.e. small centrifugal effect; at first order the flow is controlled by the 
linear equations: 



Rot (e z A u — re r 9i - EAU) 
n 2 u r = E T A@i 



(4.5) 



div u = 
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It is quite easy to see that if viscosity is neglected, the thermal wind 

u = (— s / n 2 (r)dr + F(s)) e v , 

(4.6) 

9i = 

is a solution of the equations. Here s = r sin 9 is the radial cylindrical coordinate 
and F(s) an arbitrary function describing a pure geostrophic solution. We thus 
find again that in such a baroclinic steady flow, no meridian circulation is allowed 
if viscosity is zero. 

The dependance between viscosity and meridian circulation may be made even 
more explicit if we project the system (14.51) on the spherical harmonics and have 
a look to first terms. Following Ricutord (198JJ), we set 



+ oo +1. 

^ £ uUiir+visr+wifr, 

1=0 m=-l 



©i = E E 07 

1=0 m=-l 



where 



R'i 



Y m p 



5 m 



vr; n , 



TP 



V x m 



with YI" 1 being the usual normalized spherical harmonic function. The projection 
of the equation leads to 



l—x dr \ r L 2 } t+1 or 



EA,wL 







-1-2 JL I r i+i w i+i 



+l(l + l)9 l m -6rn 2 (r)N 2 5 v . 



n 2 (r)u l r 







where N 2 = ^ y and A\_ ll B\_ x , ... are coupling coefficients. This is the set of 
ordinary differential equations which can be solved numerically to obtain the full 
solution. However, for the pedagogical matters we are interested in, it is sufficient 
to consider £ < 2. We thus have: 



EAmj 1 + A^r- 3 -^ ( r 4 u 2 ) = 



EA 2 A 2 (ru 2 ) = Bfr 1 -^-! ] + 66> 2 - 6rn 2 7V ; 



Ej* A 2 2 



n 2 u 2 = 



In the axisymmetric problem which is considered here u e 's represent the radial 
velocity and therefore meridian velocity while w 's describe the azimuthal velocity. 
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The first equation shows that if viscosity is zero, i.e. E — 0, then u 2 = 0, hence no 
meridian flow occurs and the temperature field is in equilibrium; the gradient of 

the azimuthal flow, BfrJ^ (^~T^J J us * compensates the baroclinic torque — rn 2 6N2- 

From these equations we may see that, when viscosity is taken into account, 
the shear drives a meridian circulation which in turns drives thermal imbalance ! 



4.6 Some additional remarks 

The foregoing example may be completed by the work of IGaraudl l|2002fl which 
presents numerical solutions of the same problem but in a more realistic set-up: 
there, the Boussinesq approximation is relaxed; instead the anelastic approxima- 
tion is used to filter out sound waves and permits the existence strong density 
contrasts of the background. Thus, radial profiles of pressure, temperature or 
de nsity, are more stellar or solar like. 

IGaraudl l|2002ft shows that steady flows are controlled by the non-dimensional 
ratio A = VN 2 /VL 2 = V/rj (see (|4.1|l for the definition of 77). For A > 1, i.e. for 
slow rotation, it is shown that the scaling of meridian velocity is indeed that of 
Eddington-Sweet, this circulation is confined to a layer close to the outer boundary. 
In the case of a fast rotation A <C 1 however, Garaud finds that meridian flows 
scale as v jr in agreement with our foregoing discussion. 

Could it be that Sweet's result applies to slowly rotating stars? In this case 
indeed the result found by Garaud gives an amplitude of the circulation which is 
independent of viscosity and which corresponds to a turn over time equal to the 
Eddington-Sweet time. Let us first note that extremely slow rotation are needed; 
indeed, A ^ 1 means r\ <C V, while we know that V <C 1 in a radiative zone. 
Hence, this result applies only when Eddington-Sweet time is at least V~ x greater 
than thermal diffusion time (i.e. Kelvin- Hclmoltz time); in such conditions, it is 
clear that Eddington-Sweet time is much longer than the actual lifetime of the 
star. Therefore, back to Busse's result, it is also clear that such a steady state is 
never reached. 



4. 7 Unsteady situations 

This discussion brings us to the question of the role of initial conditions. Until 
now we focused our attention to the case of steady flows, however if the star is 
initially rotating slowly so that Eddington-Sweet time exceeds nuclear time, it is 
clear that, according to Busse, the rotational state of the star will be marked by 
the initial conditions. In such case modeling the actual dynamical state of star 
requires the knowledge of an enormous number of parameters, namely the initial 
velocity field (at least!). 

However, as shown by observations, young stars are rather rapid rotators and 
therefore we may expect that they reach a sort of universal state of differential 
rotation. But it is also well known that in the course of their lives, stars lose 
angular momentum and thus progressively slow down their rotation. 
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4.8 Spin-down processes 

This remark brings us to the r ole of winds and angular momentum losses. As 
already stressed by Zahn (1992), this is certainly the main source of mixing for 
main sequence stars. The supposed universal differential rotation which comes out 
of the balance between baroclinic torque and (turbulent) viscosity is thus perturbed 
by the shear produced by the angular momentum loss. In fact this becomes the 
dominant source of turbulence in radiative zone and the dynamical state of the star 
is now controlled by the mechanisms which insure the angular momentum transfer 
between the different layers. It may be turbulent diffu sion, but also gravity waves 
emitted by a neighbouring convective zone (see iTalon & Charbonnel [2003) . 

To give an intuitive idea of how a spin down process works, it is useful to go 
back to the laboratory. Consider a fluid contained in some rotating container; at 
some time when the fluid is in perfect solid rotation, the rotation of the container 
is slowed down by some small amount (so that all the spin-down process remains 
controlled by linear equations). After some rotation, Ekman boundary layers 
form on the wall of the container which are not parallel to the rotation axis. 
These layers are the seat of shear flow parallel to the wall and when the wall 
imposes a no-slip boundary conditions, dissipation is very intense there. However, 
in addition to this dissipative process, Ekman layers also pump the fluid into or 
out of the layer. By this process, the layer drives a circulation in the bulk of the 
fluid mass, thus inducing a transport of angular momentum from the bulk to the 
layer. In a laboratory experiment the friction of the fluid, through viscosity, on 
the container's wall, finally removes the angular momentum from the fluid. The 
circulation, induced by the Ekman layer, plays a fundamental role in the spin-down 
process since it governs the time-scale of the process. This time scale is indeed 

much smaller than the diffusion time scale L 2 jv. 

Back to the case of stars there are some differences, which still need to be 
investigated in detail, but the gross features are certainly quite similar. 

There are of course no rigid wall at the surface of a star; however, if we consider 
the case of a solar type star with an outer envelope threaded by magnetic fields, it 
is likely that the convective envelope will lose quite rapidly its angular momentum 
and behave like a rigid body stressing the lower radiative zone. There everything 
depends on the viscosity : if strong enough it can extract angular momentum at the 
desired rate imposed by the wind. On the contrary, the shear layers may become 
unstable and turbulence should develop enhancing the viscosity; the foregoing 
situation may be restored and extraction takes the imposed rate; however, shear 
has a maximum value and therefore the stress (equal to shea r times visco sity) 
may not be high enough; in such a case, called strong wind in IZahnl l)l992j) . the 
radiative zone will slow down at a rate controlled by turbulent diffusive processes. 

In this description we have not described the effect of stratification and heat 
diffusion which make the whole picture a little more complex. In fact the whole 
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problem has been only touched upon and the whole non-linear scenario is waiting 
exploration . Basi c re ferences on spin -down process in the laboratory context are 



and small Prandtl number can be found in [Fricdlandcr (1976); an application to 
the stellar (solar) context has been tempted bv lZahn et alJ l|l997l) . 

5 Summarizing essential points 

To end this lecture I would like to stress the most important points, those one 
should keep in mind when dealing with rotating radiative zones. 
First we need to stress that 



is not an equation determining the velocity field. It determines the temperature 
once v is given by the dynamics; there is no other mean to determine the temper- 
ature. Giving an a-priori-approximate-solution for T and using it to derive v is 
just an incorrect reasoning. Any approximate solution of v needs to be a solution 
of both the momentum and energy equation (or a combination of both). 

It is another error to specify to rotation state (usually assumed rigid) as it 
will quickly (on a Kelvin-Hclmoltz time however) evolve to a differential rotation 
state thanks to the diffusion of heat. Even if rotation is slow, such a reasoning 
is misleading since the rotation state at time t is the result of an evolution which 
cannot be discarded; initial conditions, which are usually a fast and differential 
rotation are not forgotten by a presently slowly rotating star. 

Finally, we have shown that viscosity (or any source of diffusion of momentum) 
plays a crucial part in the dynamical evolution of a star, in particular in mixing 
processes. An isolated star (no wind) cannot sustain a large-scale circulation in a 
quasi-steady state if no diffusion process occurs. 

Finally, let us stress that the loss of angular momentum is another crucial 
aspect for the mixing processes since these losses also drive meridian circulation. 
It underlines again the importance of the dynamical history of a star when one 
tries to understand its present surface abundances. 

I would like to thank Jean-Paul Zahn for the many interesting discussions we had on this subject 
and his kind reading of this manuscript. 
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